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I propose a new mechanism of intense high-frequency electromagnetic wave generation by spatially 
uniform stacked Josephson junctions at zero magnetic field. The ac- Josephson effect converts the 
dc-bias voltage into ac-supercurrent, however, in the absence of spatial variation of Josephson phase 
difference, does not provide dc-to-ac power conversion, needed for emission of electromagnetic waves. 
Here I demonstrate that at geometrical resonance conditions, the spatial homogeneity of phase can 
be spontaneously broken by appearance of breathers (bound fluxon-antifluxon pairs), facilitating 
effective dc-to-ac power conversion. The proposed mechanism explains all major features of recently 
observed THz radiation from large area Bi2Sr2CaCu2 0s+a; mesa structures. 



Auto-oscillations are self-sustaining oscillations in non- 
linear systems with frequency and often amplitude inde- 
pendent of driving force. Oscillatory behavior in plasma 
yj, operation of a clock and multivibrators, voice and 
sound of musical instruments, all are examples of auto- 
oscillations. Typically, periodic auto-oscillations are ex- 
cited by a constant force. The frequency is determined 
either by an internal resonance (e.g. a cavity mode), 
or a characteristic relaxation time. Auto-oscillations are 
widely used for generation of microwaves, e.g. in Gunn 
diodes, magnetrons and klystrons. 

Recently a significant THz emission has been reported 
at zero magnetic field from large area E^S^CaC^Og+x 
(Bi-2212) mesa structures [I, |3j, which represent natu- 
ral stacks of atomic-scale intrinsic Josephson junctions 
(IJJs). Although the emission occurs at geometrical 
(Fiske) resonance frequencies, association with conven- 
tional Fiske steps 0| is problematic because their ampli- 
tude is zero at H — 0. Furthermore, significant emission 
power would require a large quality factor Q ^> 1 [5|. 
Since Q of Fiske steps is inversely proportional to the 
mesa size [j, it is not clear how it could be sufficiently 
large for large mesas at large operation temperatures and 
quasiparticle damping due to large self- heating. There- 
fore, the observed radiation is quite puzzling and remains 
a mater of intense discussion [5l4l2j . 

In this letter I propose a new mechanism of emission 
from a spatially uniform stack of Josephson junctions at 
zero applied field. It is shown that at geometrical reso- 
nance conditions the homogeneous state becomes unsta- 
ble with respect to formation of a breather lattice in the 
stack. Breathers are bound fluxon-antifluxon pairs [13j. 
They effectively couple the dc-bias to the ac- Josephson 
oscillations and makes resonances self-sustainable, once 
ignited. It is argued that the discovered breather auto- 
oscillations explain experimental features of the zero-field 
emission from large Bi-2212 mesas. 

Radiation from large IJJs is following the ac- Josephson 
relation Q, Q. Although the ac- Josephson effect con- 
verts the dc-bias voltage into the oscillating supercur- 



rent, this does not ensure radiation from junctions be- 
cause the supercurrent is non-dissipative and can not by 
itself pump energy from the dc-power supply into radia- 
tion. Such a dc-to-ac power conversion can be achieved 
via the Lorentz force 

F L = sIx B, (1) 

(per junction), where I is the dc-bias current and s is the 
stacking periodicity (s ~ 1.5 nm for IJJs). Therefore, 
emission by means of the ac- Josephson effect requires fi- 
nite Bi(x) 7^ 0, which is connected with the finite Joseph- 
son phase gradients ^ (i is the junction index). 

Structural inhomogeneity may lead to some coupling 
of the dc-bias to the ac-resonance field at H — [7|, [l4j . 
The inhomogeneity can be caused by variation of the crit- 
ical current density, or nonuniform bias current distri- 
bution [Isj. In large Bi-2212 mesas the inhomogeneity 
can also be induced by uneven self-heating at large bias 
[9] and by defects in Bi-2212 single crystals. Associated 
phase gradients can be attributed to self-field caused by 
nonuniform current flow. This may lead to appearance 
zero-field Fiske steps [14J. However, the observed emis- 
sion is hardly explained by zero-field Fiske steps: large 
flux quantization field in atomic scale IJJs would require 
very large inhomogeneity [7j. Furthermore, such inho- 
mogeneity could hinder mutual synchronization of IJJs 
[4J and suppress collective Fiske resonances, required for 
coherent amplification of the emission power. 

Fluxons create a large phase gradient ~ 7t/Xj (Aj 
is the Josephson penetration depth) and thus provide 
an effective dc-to-ac-power conversion. Usually fluxons 
are introduced by applying in-plane magnetic field. In 
long Josephson junctions, L ^> Aj, they can be trapped 
at if = 0. Emission scenarios, involving quasistatic 
semi-fluxons/antifluxons [6[, or fluxons/antifl uxons |8[ in 
mesas, were proposed. For a single junction static fluxon- 
antifluxon pairs are unstable, because they tend to col- 
lapse and annihilate. But for stacked junction fluxon- 
antifluxon pairs in neighbour junctions are stable [HIQ- 
Fluxon-antifluxon modes at H = has been clearly ob- 
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served in large Bi-2212 mesas, as they cause multiple val- 
ued critical current [16j . Fluxons experience the Lorentz 
force, which facilitates efficient pumping of dc-power into 
kinetic energy of fluxons. Upon collision and annihilation 
of fluxons at junction edges, radiation pulses are pro- 
duced. The corresponding flux-flow emission is well stud- 
ied for single junctions [18]. Significant flux-flow emission 
from stacked junctions requires both stabilization of the 
square fluxon lattice and high quality geometrical reso- 
nances [Hj], because the emission power is P ra d oc Q 2 and 
the linewidth oc Q~ l . For IJJs, this can only be achieved 
at strong fields H > <£> /\js ~ 2T, and for small Bi-2212 
mesas [4|, LL9j . None of those requirements is realized in 
case of zero- field emission from large mesas. 

Several non-Josephson emission mechanisms at H — 
are also known. Stable fluxon- ant ifluxon pairs in stacked 
junctions can be unbound by a large enough transport 
current and may start shuttling in the stack, leading 
to appearance of zero-field steps (ZFS) in I — V curves 
[l?}- Some emission of electromagnetic waves does oc- 
cur at ZFS 0, Ull , however, it occurs at subharmonics 
of the Josephson frequency [20| and the emission power 
is small because the fluxon is not annihilated upon col- 
lision with the edge, but is reflected as an ant ifluxon. 
Therefore, unlike in the flux-flow case, only a minor part 
of fluxons energy can be radiated. Moreover, emission 
power at ZFS can not be large because the fluxon will 
not be reflected if it will loose significant part of its en- 
ergy upon the collision with the edge [Hj. In stacks, ZFS 
can be accompanied by non-Josephson Cherenkov radia- 
tion (22I, due to partial decomposition of fluxons into 
travelling plasma waves [23], |24|. Also, recombination of 
injected nonequilibrium quasiparticles leads to genera- 
tion of bosons 
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26] . The nonequilibrium emission is 
direct (does not involve the ac- Josephson effect) and can 
provide very high efficacy [I?}. It is most effective at 
H = because it benefits from the sharp gap singularity 
in the quasiparticle density of states. 

Analysis of emission from stacked Josephson junctions 
requires implementation of proper radiative boundary 
conditions into the coupled sine- Gordon equation. In the 
numerical simulations, presented below, non-local radia- 
tive boundary conditions, derived in Ref. [5], are em- 
ployed. The emission is facilitated by the finite radiation 
impedance Z. Simulations are made for a stack of N = 10 
identical, uniform junctions with parameters typical for 
optimally doped Bi-2212 IJJs [4|, [l9j. The damping pa- 
rameter was varied from strongly underdamped a — 0.01 
to overdamped a = 1. Note that IJJs exhibit hysteresis 
in I-V characteristics, i.e., remain underdamped, almost 
up to T c [28j. The emission power is scaling with the 
width of the stack w and is normalized to w — 20/im. 

Figure 1 (a) shows I-V characteristics for stacks with 
the length L — 20Aj for a — 0.01 (black, gray and blue) 
and 0.1 (red, magenta). Different colors represent simu- 
lations with different initial conditions, corresponding to 
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FIG. 1. (Color online), (a) Normalized I-V characteristics 
for long, uniform, underdamped stack Josephson junctions at 
zero applied field. Appearance of pronounced steps, caused by 
breather auto-oscillations, is seen. Different colors correspond 
to simulations with different initial conditions with certain 
amount of trapped fluxons/ ant ifluxons. Inset shows zero-field 
steps at low voltages, (b) Normalized radiative power along 
the same I-V's. Pronounced maxima occur at collective in- 
phase cavity resonances, indicated by grid lines. 



certain amount of trapped fluxons and antifluxons. At 
low bias this leads to appearance of ZFS [17|, shown in 
the inset. At larger bias, distinct resonances appear, re- 
sembling Fiske steps. 

Figure 1 (b) shows the radiation power P ra d (from one 
edge), normalized by the radiation impedance Z. All 
presented simulations are made for very large Z, so that 
radiative losses are much smaller than internal resistive 
losses. Under those circumstances, the product P ra dZ is 
constant [5]. From comparison of Figs. 1 (a) and (b) it is 
seen that some steps are accompanied by strong emission, 
others to the contrary, correspond to emission minima. 
Distinct emission maxima occur at voltages of collective 
in-phase Fiske steps, shown by the grid in Fig. 1 (b). 

Figure 2 shows snapshots of spatial distributions of 
phases tpi (top panels), magnetic inductions Bi (middle 
panels, ~ 0.184Gs per division) and ac-components of 
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FIG. 2. (Color online). Snapshots of phase (top panels), mag- 
netic induction (middle panels) and ac-component of volt- 
age (bottom panels) at four different breather resonances 
for L — 20Xj overdamped a = 1 (a) and underdamped 
a = 0.1 (b-d) stacks. It is seen that electromagnetic field 
forms standing-wave patterns in the stack. Resonant cavity 
modes (m, n) are indicated in the figures. 



voltages Vi in each junction (bottom panels, ~ 209/iV 
per division) at for different resonances for (a) over- 
damped a — 1 and (b-d) underdamped a = 0.1 stacks 
with L — 20 Aj. The junction color code is represented 
by the dc- voltage annotation in Fig. 2 (a). It is seen 
that electromagnetic field forms two-dimensional stand- 
ing wave patterns in the stack, similar to Fiske resonances 
y, |5|, |29j. Resonant modes (m,n) are characterized by 
the wave numbers k a b — mn/L in-plane and k c — nn/Ns 
in the c-axis direction [29j. 

From Fig. 2 it is seen that there is no increment 
of phase in the junctions ipi(x = 0) = y>i(L), which 
means that the net magnetic flux in each junction is 
zero. The observed modulation of phase is, therefore, 
caused by a spontaneous formation of ordered breather 
lattice, consisting of similar breather chains in each junc- 
tion. Breather chains couple to the dc-power supply via 
the Lorentz force and effectively pump energy into ac- 
Josephson oscillations. This makes breather resonances 
self-sustaining, once ignited. The frequency of oscilla- 
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FIG. 3. (Color online). Spectra of voltage oscillations in the 
middle junction for overdamped stacks, a = 1, with L — 20Aj 
and 100Aj. Stacks are biased at the same dc- voltage, corre- 
sponding to in-phase cavity modes (2,1) and (10,1), respec- 
tively. Spectra are similar and have narrow maxima at the 
Josephson frequency cuj. Note that amplitudes and linewidths 
of the two resonances are similar, despite different L. 



tions is adopted to the nearest cavity mode in the stack. 
Such behavior is typical for auto-oscillation phenomena, 
as mentioned in the introduction. Therefore, I refer 
the discovered new type of zero-field resonances to as 
breather auto-oscillations. 

The breather amplitude is not quantized and the phase 
variation can acquire any value in the range — 2tt < 
Atpi < 2n. Fig. 2 (a) represents a special case, when the 
phase amplitude at a certain time is close to ±7r, similar 
to the case discussed in Ref. [6]. Fig. 2 (d) shows another 
special case when fluxons and antifluxons in the breather 
are well separated and the amplitude Atpi ~ ±27r, similar 
to the case considered in Ref. |8[- However, even in those 
cases, it should be realized that A(fi are not constant but 
are oscillating in time. In general the amplitude Aipi can 
be arbitrary, as shown in Figs. 2 (b,c). 

Figure 3 shows the spectra of voltage oscillations in 
the middle junction % — 5 at the edge of the stack x — 0. 
Frequency is normalized by the Josephson plasma fre- 
quency uj p . Results are shown for two overdamped stacks 
a = 1 with different lengths L = 20Aj (solid line) and 
L — 100Aj (dashed line), biased at the same relative 
current I/I c — 1.4 and having the same voltage. For 
L — 20 A j it corresponds to the in-phase resonance (2,1), 
shown in Fig. 2 (a), and for L — 100Aj to the in-phase 
mode (10,1), with the same spatial separation between 
nodes. It is seen that in both cases the radiation spectra 
are similar and consist of a sharp maximum at the pri- 
mary Josephson frequency, marked by the upward arrow, 
and several small harmonics. 

Simulations in Fig. 3 demonstrate that breather auto- 
oscillations are not hampered even in the overdamped 
case a — 1 and in very long stacks. Furthermore, neither 
the amplitude, nor the linewidth of oscillations is dete- 
riorated in the longer stack. This implies that quality 
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factors of the two resonances are similar. Quality factors 
of geometrical resonances were considered in Ref . [4j : for 
the mode (m, n), Q m ,n — rrnr(c n /co)(\/ L) /a, where c n 
is the velocity of mode n and cq is the Swihart velocity of 
a single junction. The in-phase mode n = 1 is the most 
important for achieving coherent high power emission. 
Using the approximate expression for c\ [4], we obtain 

W2(iV + l) 
(L/Xj)a 

valid for N < 7r\ a b/s ~ 400, where \ a b ~ 200nm is the 
London penetration depth. It is seen that modes with the 
same L /m from Fig. 3 indeed have the same Q~y/2> 1, 
despite junctions are long and overdamped. The quality 
factor of the junction, Qq = Q(oj — uj p ) = 1/a, should 
not be confused with the quality factor of geometrical 
resonances Q m ,n = (^m,n/^ P )Qo- 

In conclusion, a new mechanism of zero-field radiation 
from stacked Josephson junctions is proposed, via spon- 
taneous breather auto-oscillations at geometrical reso- 
nance conditions [30]. It explains all major experimental 
features of THz emission from large Bi-2212 mesas [2|, |3j: 

i) Breather auto-oscillations lead to powerful coherent 
radiation from uniform mesas at zero magnetic field, be- 
cause breather can effectively couple ac- Josephson oscil- 
lations to dc-power supply via the Lorentz force. 

ii) They excite conventional cavity modes in the stack 
with the frequency following the ac- Josephson relation 
and inversely proportional to the junction length. 

iii) Coherent flux-flow emission requires stabilization of 
the rectangular (in-phase) fluxon lattice, which is coun- 
teracted by fluxon-fluxon repulsion [ljj. Establishment 
of in-phase breather oscillations is much easier because 
of much smaller breather-breather repulsion. 

iv) Unlike conventional Fiske resonances, which have 
only one strong mode for a given field (the velocity 
matching mode with two nodes per fluxon) |4| , the wave- 
length of breather chain can be arbitrary. This flexi- 
bility provides an important advantage with respect to 
Fiske resonances, because it allows maximum amplitude 
for any geometrical resonance at H — 0. 

v) The flexibility of breather resonances ensures that 
some, high enough, in-phase resonances would have large 
quality factor, irrespective of the damping parameter or 
the junctions length. This allows intense coherent emis- 
sion of electromagnetic waves with narrow linewidth even 
in very large mesas and at elevated temperatures. 

vi) Emission due to breather auto-oscillations is sup- 
pressed by small in-plane magnetic field. Numerical sim- 
ulations (not shown) demonstrated that with increasing 
field, breather resonances are gradually substituted by 
conventional Fiske resonances, which have smaller am- 
plitude at low fields. 

vii) Breather resonances depend on initial conditions, 
which resembles metastability of emission from Bi-2212 



mesas. Essentially breather auto-oscillations need to be 
ignited to become self-sustainable. Here they were ig- 
nited with the help of trapped fluxons/antifluxons. In 
experimental situation it is likely that the ignition is facil- 
itated by the non-equilibrium mechanism [25H27| . which 
is most intense at large injection current densities just 
beneath electrodes, as seen in laser microscopy |3j. 

viii) From Eq. (2) it follows that coherent emission 
is improved in mesas with larger number of IJJs, up to 
TV ~ 400, which is the optimal number for achieving large 
emission power, without too large self-heating. 
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